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The population activity of random networks of excitatory and inhibitory leaky integrate-and-fire 
(LIF) neurons has been studied extensively. In particular, a state of asynchronous activity with low 
firing rates and low pairwise correlations emerges in sparsely connected networks. We apply linear 
response theory to evaluate the influence of detailed network structure on neuron dynamics. It turns 
out that pairwise correlations induced by direct and indirect network connections can be related to 
the matrix of direct linear interactions. Furthermore, we study the influence of characteristics of the 
neuron model. Interpreting the reset as self-inhibition we examine its influence, via the spectrum 
of single neuron activity, on network autocorrelation functions and the overall correlation level. 
The neuron model also affects the form of interaction kernels and consequently the time-dependent 
correlation functions. We finally find that a linear instability of networks with Erdos-Renyi topology 
coincides with a global transition to a highly correlated network state. Our work shows that recurrent 
interactions have a profound impact on spike train statistics and provides tools to study effects of 
specific network topologies. 



I. INTRODUCTION 

The dynamics of networks of interacting units, e.g. 
correlation and synchronization, depends on the prop- 
erties of the nodes, as well as on the organization of 
their interactions. Neurons in cortical networks typically 
spike asynchronously despite a large number of synapses 
coupling the neurons together. This behavior has been 
explained by random balanced networks of excitatory 
and inhibitory neurons [TJ [2]. Later studies examined 
in greater detail how connections between neurons me- 
diate correlations [3], but inhibitory network feedback 
decorrelates activity [I] . Residual correlations are shaped 
by the detailed structure of the network and potentially 
influence neuronal information processing. Because the 
spiking operation typically incorporates a thresholding 
mechanism, single neuron models are intrinsically non- 
linear. A further complication results from the fact that 
even in sparse networks recurrent interactions or feed- 
back loops have a large influence on dynamics, leading to 
chaotic dynamics [TJ[5]. This implies that analytic results 
on a network level are difficult to obtain. Here we use 
the impulse response to construct an efficient lineariza- 
tion in the asynchronous-irregular state of a network of 
LIF neurons. Linear models represent a powerful tool to 
study structured networks in detail [BJ [7j- As it turns 
out, effects of the network can be well described in such 
a linear framework. This provides a convenient way to 
include the full recurrent connectivity of the network and 
examine its interaction with single-neuron dynamics and 
demonstrates that the implications of recurrent feedback 
can be understood as a linear, albeit not trivial effect. 
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II. NETWORK MODEL 

The LIF neuron is a standard model in studies of neu- 
ral network dynamics. The dynamics of the membrane 
potential V k of neuron k is governed by the equation 

T m Vk = -V k - T m VgS k + T m JkjSj, (1) 

j 

where r m denotes the membrane time constant and J^j 
the synaptic weight between presynaptic neuron j and 
postsynaptic neuron k. A spike is emitted when the 
membrane potential exceeds a threshold Vt- Spikes are 
modeled as <5-shaped events, such that the spike-train of 
neuron j is Sj(t) — S(t — tj). The membrane potential 
is decreased by an amount Vq after each spike. We sim- 
ulate random networks comprising N neurons of which 
e = 80 % are excitatory and 1 — e = 20 % are inhibitory. 
Connections are established independently with proba- 
bility p = 0.1, such that a diluted network is obtained. 
Synaptic weights are Je for excitatory and Jj = —^Je 
for inhibitory neurons. A value of 7 > 4 then implies 
a relative dominance of inhibition in the network. To 
study the system in a regime where neurons are mainly 
driven by membrane potential fluctuations, we assume 
that all neurons in the network receive Poisson-like in- 
put from external neurons with synaptic weights Je and 
rate y ex t — Vr/(JET m ), such that the mean membrane 
potential is just high enough to reach threshold. In sim- 
ulations we use T m = 20 ms, a refractory time r r = 2 ms 
after each spike, a threshold Vt = 20 mV and a reset by 
Vg = 20 mV if not mentioned otherwise. Network simu- 
lations were conducted with the NEST-simulator [5]. 

The dynamic state of random networks of LIF neurons 
depends, besides the strength of the external input, on 
the ratio of excitation and inhibition. We consider neu- 
rons in balanced or inhibition dominated networks which 
spike asynchronously and irregularly at stationary firing 
rates, see Fig.HTa). As a measure of correlations between 
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FIG. 1. (color online) (a) Spike trains in a random LIF 
network (N — 1250,7 — 5, Je = 0.1 mV, simulation time 
5 • 10 4 s). Covariances and fluctuations of population spike 
counts (bottom) are small, (b) Autocovariance of neurons 
with low and high rate. Normalization over bin size At was 
chosen such that the integral of the peak corresponds to firing 
rate, (c) Histogram of the response to external excitatory and 
inhibitory input spikes, averaged across neurons, as well as fits 
to a function cos(z/r)e _T//a (dashed lines). 



spike trains we use the covariance density function 
C kj {T) = (s k (t)s 3 (t + T)) - (s k {t)}( Sj (t)}. 



(2) 



The autocovariances C kk exhibit troughs around r = 
because the reset after each spike leads to an effective 
dead time, panel (b). The effect of a single input spike 
of weight J on the firing rate of neuron fc, the impulse 
response gk(J), panel (c), can be measured as the nor- 
malized covariance of some external Poisson input s e (t) 
coupled to k with weight J, 



(Se(t + T),S k (t))-( Se (t))(s k (t)) 
(Se(t)) 



(3) 



III. LINEARLY INTERACTING POINT 
PROCESSES 



where Cq(t) = C(t) — YS(t) denotes the continuous part 
of the covariance matrix. The rate matrix Y is diagonal, 
with elements \Y]kj — Vk^kj, coinciding with the inte- 
grated autocovariances of Poisson processes. This rep- 
resents the fact that isolated nodes with no interactions 
are modeled as Poisson processes. In the Fourier domain, 
the solution of ^ is 



C(w) = [1 - G(oj)]~ 1 Y[1 - G*(lj)]~ 



(7) 



where G* denotes the transpose and complex conjugate 
of G, see [9] for details. The reset after each spike af- 
fects the autocovariance functions of LIF neurons. It is 
therefore not feasible to assume that isolated nodes spike 
as Poisson processes. Nonetheless, their activity can be 
described in the Hawkes framework, if reset effects are 
taken into account. This can be done by interpreting the 
reset as self-inhibition. 

In a network with synaptic coupling matrix G syn , self- 
inhibition is described by diagonal elements —9k, so that 
the full coupling matrix is G = G syn — O, with [Q]kj = 
SkjO k - In this case, Eq. can be rewritten as 



C(w) = [1 - G(w)] -1 Y(u;)[l - G*(u)Y 



(8) 



with G(w) = (1 + e(cj))- 1 G syn (o;) and Y(u) = (1 + 
0(w)) _1 F(l + 0*(w)) _1 . One can interpret G as effec- 
tive interaction matrix between neurons with reset. The 
elements of Y coincide with the autocovariances of uni- 
variate Hawkes processes with self-inhibition, which are 
more regular than pure Poisson processes. Covariances 
decrease, because the couplings [G syn ]kj are weakened to 
[G syn ]kj/ (1 + 9 k ). 



IV. MAPPING TO LIF NETWORKS 



A simple description of spike train dynamics is pro- 
vided by the model of interacting point processes intro- 
duced by Hawkes [S]. Spike trains Sk(t) are interpreted 
as realizations of Poisson processes with dynamic rates 
Uk(t) = ( s k{t))- The spike train vector s(t) and the rate 
vector y(t) are used to describe a population of neurons. 
External input is assumed to evoke a baseline firing rate 
2/0- Presynaptic spikes at time t — r cause a transient 
response. All interaction kernels are collected in the ma- 
trix G(t), with G(t) = for r < 0. The dynamics of the 
system is then defined by the equation 

/oo 
G(t-T) S (T)dT = y Q + (G*s)(t). (4) 
-OO 

In a stationary system, the time averaged rates are 

y=( 1 -J G(r)dT^j y , (5) 

where y = (s(t)) was used. The matrix of crosscovari- 
ances C (r) is determined by the self-consistency equation 

Co(t) = G(t)Y + (G*C )(t), (6) 



The prediction from the linear theory in ([8]) has been 
rewritten so that all quantities are experimentally acces- 
sible and can therefore be tested in simulations. In a 
first step, we show how the linear response of neurons 
in a network can be identified with the elements of the 
effective coupling matrix G. The coupling matrix deter- 
mines the extent and distribution of covariances in the 
network. In the second step, the effective autocovariance 
can be determined from measured autocovariances for a 
quantitative prediction of crosscovariance functions. 

We shortly summarize some definitions at this point. 
The effective coupling matrix G results from reset effects 
in the actual coupling matrix G. External input with 
coupling g ex t(u>) to linear networks or synaptic weights 
Je,Ji to networks of LIF neurons can be used to in- 
fer an estimate g(ui) of the impulse response of network 
neurons. Integrated estimated impulse responses will be 
denoted by gE,gi for excitatory and inhibitory neurons, 
respectively. 

Analogously as in ([3]), the impulse response gk{gext) 
of neuron k is measured as the covariance between the 
spikes of an external input node coupled to k with kernel 
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g ex t- This covariance can be calculated using Equation 
([8]), where the coupling matrix has to be extended to 
include the external input. Specifically, this is done by 
calculating covariances in a network with coupling ma- 
trix G' with entries G' kl = g ext , = Gkj for 
k,j = l...N and otherwise, such that the new node 1 
corresponds to the external input. The result is 



gk(g e xt) 



2/H 



gext 

(i + o k ) 



{(l-G)-% k . 



(9) 



Here, j/n is the rate of the external input. The factor 
[(1 — G) _1 ]fcfc relates to the fact that, as the neuron is 
embedded in a network, the total response that is mea- 
sured includes the direct response as well as the indi- 
rect response via its postsynaptic partners, which in turn 
project back to the original neuron. For large diluted 
networks with (Gkj) <S 1 this factor is negligible. There- 
fore, gk(gext) ~ 9ext/0- + 8k)- If gext is chosen identical 
to the internal couplings, measurement of the response 
then returns the elements of the effective coupling ma- 
trix. In simulated networks of LIF neurons, measured 
impulse responses can thus be identified with effective 
couplings, [G]kj — gk{Jkj)- As in our simulations synap- 
tic weights are either Jg or J/, the responses (^(Je) and 
gk(Ji) have to be measured. If the synaptic connectivity 
Jkj is known, the coupling matrix is then fully deter- 
mined. If only certain statistics of the covariances are of 
interest, knowledge of statistical properties of the con- 
nectivity matrix may be sufficient, see Section |VIII| In 
principle, gk{J) is different for each neuron fc, as vari- 
ances in input affect the dynamical regime it is in, and 
has to be determined individually. In random networks, 
the variance in input is relatively small, and to reduce 
simulation costs, we assume identical impulse responses 
for all neurons. The response can then be measured as 
an average across a small population of probed neurons. 

Because the effective autocovariances Y are not di- 
rectly observable, the autocovariance functions of the 
nodes, Ckk have to be measured instead. According to 
0, they are linked to the Y by C kk (ui) = £\ |[(1 - 

G{ui))^]\ 2 yjj{Lo). Inverting this equation removes the 
network feedback. Together with the coupling matrix, 
this allows for a linear prediction of the remaining cross- 
covariances. 



V. COVARIANCES IN RANDOM NETWORKS 

The linear model then predicts covariance functions 
between individual neurons in frequency or, equivalently, 
in time domain, see Fig. [2j Predictions are particularly 
accurate for lower frequencies. 

A reset after each spike is a feature of the biophysi- 
cal mechanism of spike generation, and an integral part 
of many neuron models. By varying the parameter Vg, 
the effective resets k are changed and its consequences 
for the covariances in the network can be studied. In 



particular, a divergence between (y) and (y) is seen for 
low values of Vg, that is small resets. This illustrates 
the point that for the estimation of crosscovariances, it 
is essential to remove the feedback from the autocovari- 
ances, and use the effective values instead. Simulation 
results are in line with expectations for the dependence 
on self-inhibition from the linear theory, Fig. [3] An in- 
crease of Ve causes weaker impulse responses and smaller 
autocovariances as well as weaker covariances. 



VI. LINEAR STABILITY 

A necessary condition for the stability of the linear 
system Q is that the real part of all eigenvalues of the 
matrix jG(t) dr is smaller than 1. This condition is also 
sufficient if the coupling kernels have the shape of a de- 
caying exponential [7]. We examined how the spectrum 
of the integrated coupling matrix scales with the size N of 
a random network, assuming a uniform connection prob- 
ability p. 

The entries of the integrated coupling matrix can be 
computed analytically [H [TU] . Starting point is the self- 
consistent rate of the recurrent network y sc , which is 
given by the solution of the self-consistency equation 



i AV T -)i)/o a 

-=T r + e" (1 + erf(u)) du, (10) 

V J(VT-V g ~fj,)/a- 



where (i = r m (J E y, 



- pNeJEy + pN(l — e)Jiy) and 
a = T m {y ext J% + pNeyJ 2 E + pN(l - e)yJf) are the 
mean and the variance of the total input to each neu- 
ron. For small synaptic weights Je,Ji, the integrated 
impulse response to a delta-shaped current is c/e,i = 



(a) o.io 
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FIG. 2. (color online) (a) Scatter plot of covariances be- 
tween neuron pairs. Red (dark gray): Cfej(w = 0kHz), 
black: Re[Ciy](cd = 0.31kHz), cyan (light gray): lm[Ckj]{oJ = 
0.31kHz). Inset: Standard deviation between prediction and 
simulation, normalized by the standard deviation of the co- 
variance distribution, increases with frequency, (b) Typical 
examples for simulated (gray) and predicted (black) pairwise 
covariance functions in the time domain. Parameters as in 
Fig. 
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FIG. 3. (color online) Larger reset decreases (a) effective and 
network autocovariances, (b) integrated impulse response and 
(c) mean and standard deviation of distribution of integrated 
covariances in a random network. Parameters as in Fig. [T] 



/ g(J E ,i,t)dt = TmJEjjj;]^^), with 



dp a 



[1 + erf ((1/ T - n)/tr) 



(11) 



z (V T -Ve-»)/** ^ + erf yy T _y _ ^ _ 



The spectrum of the adjacency matrix of random net- 
works comprising excitatory and inhibitory neurons has 
been described previously [IT] . In the limit of large net- 
works, one eigenvalue corresponds to the average input 
to each neuron, m ~ pN(ec/E + (1 — £)<?/)• The remain- 
ing eigenvalues are located within a circle in the complex 
plane of radius p — \/N(eg E + (1 — e)gj)p(l — p). If the 
network is balanced or dominated by inhibition, m < 
and stability depends only on p. 

Fig. [4] shows that for large system size §ej decreases 
faster than iV -1 / 2 in a balanced network. As p oc A 1 / 2 , 
this implies small p and guarantees stability and low co- 
variances, in line with previous work [2]. In the case 
of inhibition dominance, the increase of input in larger 
networks depresses firing rates, causing an even more 
rapid decay of impulse responses. However, in interme- 
diate sized networks with weak inhibition and relatively 
strong synapses, the system can become linearly unstable 
as p grows above 1. The analytical values for impulse re- 
sponses and corresponding p are in good agreement with 
simulations for small networks far away from the unsta- 
ble regime. Activity in networks in the linearly unstable 
regime is bounded, as firing rates saturate e.g. due to 
the refractory time, but an asynchronous irregular state 
turns into a state with strong pairwise covariances, Fig. 
[5]ja). The average integrated covariances reflect this phe- 
nomenon, they increase steeply for large p, panel (b). 



VII. INFLUENCE OF COUPLING KERNELS 

The shape of the coupling kernels is not well described 
by a decaying exponential, compare Fig. [TJ We therefore 
take a closer look at the stability of the system described 
by Eq. Q. It can most easily be analyzed, when the 
spike train s(t) is replaced by the conditional intensity 



y(t). The resulting equation is 

y(t)=y + G*y(t). 



(12) 



A result from Paley- Wiener theory [T^] is that a linear 
system is stable if and only if the real part of all poles 
of the transfer function, which is defined in the Laplace 
domain as 



C[T](s) 



£[yo](s) 



(13) 



and relates input yo and output y, is smaller than 0. 
After a Laplace transform of Eq. ( 12 ) one obtains 



C[T](s) = (l-C[G}(s)Y 



(14) 



the typical expression for a system with feedback. An 
alternative formulation for the stability condition is that 
1 must not be part of the spectrum of C[G](s) for any 
s from the right half-plane. This condition on the spec- 
trum can be readily evaluated, when, for simplicity, we 
assume that all coupling kernels have the same time de- 
pendence, so that excitatory and inhibitory couplings 
only differ by a scaling factor. The coupling matrix is 
then of the form G — g(t)K, where K is a scalar ma- 
trix. The bulk spectrum of K is circular and has a 
radius px- The spectrum of £[G](s) then is a scaled 
version of the spectrum of K, with radius p/f|£[g](s)|, 
compare Fig. p[c), as the rotation caused by the com- 
plex factor C[g](s) does not change the shape of the bulk 
spectrum. If liniRo^^oo |£[g](s)| = 0, the bulk spectrum 
is arbitrarily small, therefore, because 1 must not to be 
part of it, stability requires that pK\C[g\{s)\ < 1 for all 
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FIG. 4. (color online) (a) Semi-analytical qe,i in random 
networks of size N for different values of Je = 0.1,0.4,0.6 
(dot-dashed, dashed, full lines), 7 = 4. (b) Bulk spectrum 
radius p of corresponding coupling matrices for 7 = 4, 4.5. (c) 
Semi-analytical couplings (dot-dashed lines) and values from 
simulations with Je = 0.1. Circles: 7 = 4.5, squares: 7 = 6. 
Empty symbols: excitatory couplings qe for 7 = 4.5, 6. Filled 
symbols: corresponding inhibitory couplings \gi\. (d) Same 
for p. 



5 




FIG. 5. (color online) (a) Average of covariance functions 
across pairs (N = 25000, Je = 0.4). (b) Integrated covari- 
ances normalized by average rates over p(N), for Je as in 
Fig. [4] Lines for different Je partially overlap. Simulation 
time 10 2 s, other parameters as in Fig. [I] (c) Eigenvalues 
of integrated coupling matrix (green/light gray) and matrix 
rescaled by g(so)/g(0) at resonance so determined from fit 
to excitatory kernel in Fig. [I] (blue/dark gray), (d) Normal- 
ized Laplace transformed impulse response |£[g(s)]|/g(0) for 
kernel fit cos(^r)e _T ' a (left) and exponential kernel (v = 0). 



s in the right half-plane. If \g{s)\ has its maximum at 
s = 0, as is the case for exponential kernels, the spec- 
trum of the integrated coupling matrix is decisive, be- 
cause £[<7](0) = J gif) dt. However, the impulse response 
of LIF neurons has resonances [13] . We approximate the 
shape of the impulse response obtained in simulations by 
the function g(r) = ae~ T / a cos(^t), compare Fig.JIJ Res- 
onances appear as maxima in C[g](s) for s ^ 0, see also 
Fig. J5jd) , left panel. These networks are therefore less 
stable than indicated by the integrated impulse response 
only. This is not the case for exponential kernels (right 
panel) . 



VIII. DISTANCE DEPENDENT 
CONNECTIVITY 

In [7] we demonstrated that, by expanding Eq. ^ into 
C( U ) = ^G"(ffl(G*rH, (15) 

covariances can be interpreted as a sum of contributions 
from direct and indirect interactions, under the condition 
that the eigenvalues of G have an absolute value smaller 
than one. In this series, the terms with {n,m} = {1,0} 
describe covariances caused by direct interactions be- 
tween nodes. The term G(ui)Y(oj)G* (cu) counts the con- 
tribution from common input and terms with higher pow- 
ers describe indirect interactions via intermediate nodes. 
To demonstrate how covariances are connected to these 
indirect interactions as well as the usefulness of the lin- 



ear approximation in networks with non-random struc- 
ture, we consider nodes arranged in a ring, where connec- 
tion probability between two neurons depends on their 
geodesic distance d. We use a boxcar connectivity func- 
tion, p a (d) = A a Q(d—a a ), where A a and a a are constants 
for a € {E,I}, meaning that the connectivity profile is 
different for excitatory and inhibitory neurons. 

In these networks the distance dependent covariance 
can be denned as c u {d) = {Ck,k+d{^))k- The quan- 
tity c u (d) depends on the average coupling 7(d) = 
tg(JE,u)pE(d) + (1 — e)g(Ji,co)pi(d) and the average 
common input 77(d) = {tg{JE,^>)) 2 PE * PE{d) + ((1 — 
t)g{Ji,uj)) 2 pi-kpi{d), where * denotes a discrete convolu- 
tion. These expressions result from the positive and neg- 
ative interactions of the two neuron populations. Anal- 
ogously to Eq. (15), c w (g0 then is decomposed in a sum 



over interactions of different orders, 



,{d) 



(d), 



(16) 



where the matrix products in ( 15 ) are replaced by convo- 



lutions of the average interactions, ci 
r y-*.^*(™- 1 ) f or 771,71 > and c^ ' 1 ^ = cL ' 1 ^ = {y)"f. Ap- 
plication of the spatial Fourier transform 



JV-l 

£ 

d=Q 



c u {d) exp(-i2irld/N) 



leads to the closed expression 



F[cJ\ = (y)(l + 



For details of the calculation compare [7]. The effec- 
tive autocovariance (y) does not affect the shape of the 
distance dependence and is merely a proportionality fac- 
tor. Application to the simulated LIF networks in Fig. 
[6] shows that the distance dependence of the covariance 
results from the sum of the direct and many indirect con- 
tributions. Both the shape and the spatial extent depend 
on indirect paths. Higher order terms corresponding to 
indirect contributions contribute strongly, compared to 
the direct interactions. For covariances on long time 
scales, represented by the zero-frequency covariance, pos- 
itive and negative contributions partly cancel, so that 
effective covariances are of moderate size. The relative 
contributions of the higher order terms vary with fre- 
quency, due to the frequency dependence of the linear 
response. As a consequence, the spatial profile of the 
distance dependence is strongly frequency dependent. 



IX. DISCUSSION 

We have shown that a linear theory describes covari- 
ances in networks of LIF neurons in approximately asyn- 
chronous states. While linearization is a common ap- 
proach to approximate complex dynamics, the existence 
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FIG. 6. (color online) Average of covariances over pairs of 
equal distance for to = kHz (a) and ui — 0.11 kHz (real part) 
(b). Covariances in the direct simulation are well reproduced 
by the linear approximation. Covariances in the linear model 
fluctuate around the analytical profile due to the random re- 
alization of the network. Contributions of different orders can 
be delineated in the analytic expression ( 16 1 (c),(d). Distance 



and frequency dependence of covariance (real part) in direct 
simulations (e) and the linear approximation (f). Simulation 
time was 2 ■ 10 4 s, p = 0.05, g = 6.5, a E = 312, 07 = 156, N = 
2500. 



sponse and interact with the network structure to de- 
fine the overall covariance configuration. In the special 
case of LIF neurons, impulse responses have been studied 
in detail [T31 HH] ■ Generally, the impulse response in a 
certain dynamical state can be measured from external 
input in simulations of more complicated neuron models 
or in experiments. Similar to the way that the fluctu- 
ation driven regime seems to favor linear dynamics for 
LIF neurons, we expect that this holds for networks of 
other neuron models as well, as long as correlations and 
population fluctuations are small. 



Although our approach is restricted to stationary con- 
ditions, it represents a first step towards a better under- 
standing of the influence of the full network structure 
on the behavior of neuron populations. One such as- 
pect that is difficult to investigate in mean field models 
is the stability constraint imposed by the radius of the 
bulk spectrum of the network related to the topological 
randomness of the system [TT]. It implies that global 
balance of excitation and inhibition is a necessary, but 
not sufficient condition for stable activity in the network. 
Simulations in the regime where stability is violated in- 
dicate that covariances become large and strong activity 
fluctuations occur. This issue has not been addressed 
much in the study of neural networks, but it was shown 
that stability constraints on growing networks can lead 
to a complex micro-structure |19j . Resonances in the im- 
pulse responses can equally affect stability, similar as in 
[2U] where bifurcations related to the coupling between 
populations were studied. 



of multiple feedback loops in recurrent networks casts 
doubts on the applicability of a linear theory. In retro- 
spect, one apparent reason is that contributions of dif- 
ferent order effectively cancel due to the interaction of 
excitatory and inhibitory populations. 

The price to pay to apply linear response theory in this 
framework is that both the linear response and the auto- 
covariances of the neurons need to be inferred from sim- 
ulations, measurements or an extended analytical frame- 
work. Even though, information about the dynamics, 
as the linear stability of the network, the distance de- 
pendence of covariances, or the relative contributions of 
indirect interactions can be inferred from the impulse re- 
sponse and the network structure alone. An additional 
asset is the possibility to compare different connectivity 
structures in a manner that does not depend on the spe- 
cific neuron model used. So far, this problem has been 
approached mainly in simulations of the full dynamics of 
the network [TH [TS] , but see [TS] , where a rate model has 
been used and |17| for an alternative linearized approach. 

Single neuron properties determine the impulse re- 



We use networks with constant connection probabil- 
ity which are diluted rather than sparse in the large size 
limit. For network sizes accessible in simulations, this 
scaling covers measured values for the estimated num- 
ber of a few thousand synapses for each neuron in local 
cortical networks where the connection probability has 
been estimated to be about 0.1. We have shown in [7] 
that the relative impact of more indirect contributions is 
more significant if the spectral radius of the network is 
large. In the present work we extend these results and 
show that, as the effective coupling strength depends on 
the frequency, the resulting differences in relative weight- 
ing of indirect contributions leads to a complex frequency 
dependence of spatiotemporal correlations. In j!6) it was 
demonstrated that small changes in the out-degree vari- 
ance can lead to a synchronous state when the system 
is close to a bifurcation, indicating that connectivity has 
a large influence if networks are almost unstable. On a 
functional level, a network where the largest eigenvalue 
of the adjacency matrix is equal to one has an optimal 
dynamic range |21j , and strong recurrence has been pro- 
posed to be a useful condition for cortical computations 

m 
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